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Abstract. Averaged inhomogeneous cosmologies lie at the forefront of interest, since 

cosmological parameters such as the rate of expansion or the mass density are to 

be considered as volume-averaged quantities and only these can be compared with 

observations. For this reason the relevant parameters are intrinsically scale-dependent 

and one wishes to control this dependence without restricting the cosmological model 

by unphysical assumptions. In the latter respect we contrast our way to approach the 

averaging problem in relativistic cosmology with shortcomings of averaged Newtonian 

models. Explicitly, we investigate the scale-dependence of Eulerian volume averages 

of scalar functions on Ricmannian three-manifolds. We propose a complementary 

view of a Lagrangian smoothing of (tensorial) variables as opposed to their Eulerian 

averaging on spatial domains. This programme is realized with the help of a global 

Ricci deformation flow for the metric. We explain rigorously the origin of the Ricci 

flow which, on heuristic grounds, has already been suggested as a possible candidate 

for smoothing the initial data set for cosmological spacetimes. The smoothing of 

geometry implies a renormalization of averaged spatial variables. We discuss the results 

in terms of effective cosmological parameters that would be assigned to the smoothed 

cosmological spacetime. In particular, we find that the on the smoothed spatial domain 

B evaluated cosmological parameters obey f% + f% + Q-g + f% = 1 , where Q-g , fig 

_A . . ~ — Q 

and fig correspond to the standard Friedmannian parameters, while f2g is a remnant 

of cosmic variance of expansion and shear fluctuations on the averaging domain. All 

these parameters are 'dressed' after smoothing-out the geometrical fluctuations, and 

we give the relations of the 'dressed' to the 'bare' parameters. While the former 

provide the framework of interpreting observations with a "Friedmannian bias" , the 

latter determine the actual cosmological model. 



PACS numbers: 04.20 , 98.80 , 02.40 K 



Regional Averaging and Scaling in Relativistic Cosmology 



2 



Introduction 

Research on cosmological spacetimes has been in the realm of general relativity for a long 
time, establishing the standard cosmological models that are based on homogeneous (and 
mostly isotropic) solutions of Einstein's laws of gravity for a continuous fluid. Spatially 
homogeneous spacetimes are understood to a high degree, and cosmologies based upon 
them certainly lie in a well-charted terrain. The difficulties or better challenges arise, 
if we want to respect the actually present inhomogeneities in the Universe. Newtonian 
continuum mechanics appears to be a simpler theory to model the inhomogeneous 
Universe, usually restricted to the matter dominated epoch on subhorizon scales. Indeed, 
most contemporary efforts for the modelling of inhomogeneities are based on Newtonian 
cosmological models. There is, however, a drawback: Newtonian continuum mechanics 
of a self-gravitating fluid is not a proper theory per se [pOfl , but has to be setup 
with suitable boundary conditions; for the cosmological modelling it is practically 
restricted to setting up the physical variables relative to a homogeneous background, 
while the (inhomogeneous) deviations thereof have to be subjected to periodic boundary 
conditions. Even though we may accept periodic boundary conditions as a necessary 
cornerstone of a Newtonian model - hence, we view the Universe as a caleidoscope 
of ever-repeating self-similar boxes that are supposed to supply a 'fair sample' - the 
introduction of a global reference background is essential to do so. This may be 
illustrated as follows. Consider Poisson's equation for the Newtonian gravitational 
potential. This potential cannot be periodic as a whole (e.g. for a homogeneous 
background it grows quadratically with distance from an origin). Moreover, solutions of 
Poisson's equation are only unique, if the spatial average of the source vanishes. Both 
requirements, periodicity and uniqueness, can be accomplished only for fields defined 
as inhomogeneous deviations from a given reference background, e.g. the standard 
FLRW models (for details compare Note that most currently employed models 

including numerical N-body simulations rest on these assumptions. These "forcing 
conditions" must be considered a drawback for the following reason: we may consider 
the spatially averaged variables as replacing the former homogeneous variables, e.g., 
the volume-averaged rate of expansion measuring the Hubble law on a given averaging 
scale. This (effective) expansion rate does not obey the Friedmann equations of the 
standard cosmological models, but the true equation features an additional source 



term due to kinematical fluctuations ||. As was first pointed out by Ellis p2| , this 
so-called "backreaction effect" is a result of the nonlinearity of the basic system of 
equations, if general relativity or Newtonian theory, lying at the heart of the problem 
how we could compare and match the FLRW standard model of cosmology with an 
averaged inhomogeneous model |E3|. It is here, were the restriction of using a Newtonian 



model becomes evident: this extra source term is a full divergence of a vector field 
and, hence, consistently vanishes on the periodic boundary. It is, however, not a full 
(three-dimensional) divergence on non-Euclidean spaces (see || for a discussion of 
this point). Hence, although it is commonly agreed that observables like Hubble's 
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"constant" or the mass density depend on the surveyed volume of space and must be 
intrinsically scale-dependent, the Newtonian models have to introduce a "largest scale" 
where these observables assume a constant value that is determined by the homogeneous 
standard cosmology and, consequently, by initial conditions given on the largest scale 
only. By construction^ an averaged Newtonian model features the characteristics of 
the standard homogeneous models. According to what has been said above, this scale 
appears artificial and we may not have such a scale where the averaged inhomogeneous 
model can be identified with or approximated for all times by the standard model. In 
other words, we may not find a global frame comoving with a standard Hubble flow and 
at the same time providing the evolution on average. As an example we point out that 
the Newtonian curvature parameter is determined by the initial data on the periodicity 
scale (e.g., a flat Einstein-de Sitter cosmology remains so during the evolution), while 
in general relativity the averaged scalar curvature is coupled to the backreaction of the 
inhomogeneities ||. Since the dynamical evolution of the curvature parameter on scales 
smaller than the periodicity scale strongly depends on the inhomogeneities (see [ lOj for 
a quantitative investigation), we can expect that a generic averaged cosmology will not 
keep the global average curvature at this initial value; it will, like any other variable, 
change in the coarse of evolution (for further discussion see ||). It should be remarked 
here that very often the argument is advanced that the backreaction term is negligible, 
because it is numerically small. On sufficiently large scales the latter is supported by 
some of the following investigations: f2|], @, §, ||, [[§], @, |39|, 0. Still, a small 
perturbation can (and as shown in |10|) will drive the dynamical system for the averaged 
fields into another "basin of attraction" implying drastic changes of the volume-averaged 
cosmological parameters although the backreaction term is numerically small. 

This "kinematical backreaction" representing, roughly speaking, the influence 
of fluctuations in the matter fields on the effective (spatially averaged) dynamical 
properties of a spatial region in the Universe, does not comprise the whole story. 
Even if we take the influence of fluctuations on the averaged variables into account, 
these variables themselves still depend on the bumpy geometry of the inhomogeneous 
averaging region. It turns out that this problem is quite subtle and lies at the heart of 
any interpretation of observables in terms of a cosmological model: observed average 
characteristics of a surveyed region are, by lack of better standards, taken as averages 
on a Euclidean (or constant curvature) space section. Any matter averaging program 
in relativistic cosmology is not complete unless we also devise a way to interprete the 
averages on an averaged geometry. The latter, however, is a tensorial entity for which 
unique procedures of averaging are not at hand. In the present paper we especially 
address this problem and propose a Lagrangian smoothing of tensorial variables as 

| Cosmologists were employing this construction for a long time without justification, and they have 
been "lucky" that the backreaction term is indeed a full divergence which in turn implies that it 
vanishes on the periodic simulation box comoving with a standard Hubble flow. Without this (non- 
trivial) property of the backreaction term cosmological N-body simulations would just be artificial 
constructions. 
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opposed to their Eulerian averaging. The present investigation will reveal a further 
shortcoming of Newtonian cosmology: curvature fluctuations turn out to be crucial and 
may even outperform the effect of kinematical fluctuations quantitatively. 

In summary we can state the following headline of our investigation of the averaging 
problem: since averaged scalar characteristics form an important set of parameters that 
respectively constrain or are determined by observations, it is a highly relevant task to 
develop a theoretical framework for averaging and scaling within which the currently 
collected datasets can be analyzed reliably and free of unphysical model assumptions. 
Newtonian models are, due to their very architecture, not free of such assumptions. In 
light of this the modelling of cosmologies has to be lifted back on the stage of general 
relativity, leaving behind the Newtonian "toy universe models", which were helpful to 
understand basic properties of the formation of structure, but have reached a dead end 
where the Euclidean periodic box is taken for real and every observational data "fitted" 
to its parameters. Fortunately, the structure of the basic equations that govern the 
averages of observables in general relativity is so close to their Newtonian counterparts, 
that it is evident to better work in the relativistic framework (compare the formal 
equivalence of the effective expansion law in Newtonian cosmology with that in 
general relativity ||). The challenge that we shall meet here, and the answer that we 
shall provide, concerns the interpretation of the average characteristics in a regional 
survey within a smoothed-out cosmological spacetime. 

Let us now come to the content of this article and describe our approach before 
we formalize it. Consider a three-dimensional manifold equipped with a Riemannian 
metric (E, <7 a &). On such a hypersurface we may select a simply-connected spatial 
region and evaluate certain average properties of the physical (scalar) variables on that 
domain such as, e.g., the volume-averaged density field, or the volume-averaged scalar 
curvature. These (covariant) average values are functionally dependent on position and 
the geometry of the chosen domain of averaging. Let us now be more specific and relate 
the averages to scaling properties of the physical variables. It is then natural to identify 
the domain with a geodesic ball centred on a given position, and introduce the spatial 
scale via the geodesic radius of the ball. We may consider the variation of this radius 
and so explore the scaling properties of average characteristics on the whole manifold. 
It turns out that this scaling depends on the intrinsic geometry of the hypersurface, 
and we shall explicitly evaluate this dependence in the neighborhood of the domain of 
averaging. Averaging regionally at every position chosen in the spatial slice at a fixed 
averaging scale, we arrive at averaged fields that, upon changing the scale, depend on 
the accidents of the regional geometry. 

We may call this point of view Eulerian - and everybody would first think of this 
point of view - in the sense that the spatial manifold is explored passively by blowing up 
the geodesic balls covering larger and larger volumes with a larger amount of material 
mass. On the contrary, the key idea of our approach consists in demonstrating that a 
corresponding active averaging procedure can be devised which is Lagrangian: we hold 
the ball at a fixed (Lagrangian) radius, and deform the dynamical variables (actively) 
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inside the balls such that the deformation corresponds to the smoothing of the fields. 
The (first) variation with scale of, e.g., the density field or the metric is then mirrored 
by one-parameter families of successively deformed density fields and metrics. We shall 
show that this deformation corresponds to a first variation of the metric along the Ricci 
tensor, known as the Ricci flow. Since this flow has received great attention in the 
mathematical literature - major contributions are due to Hamilton P7fl , [ ffijfl - we shall 
so translate the averaging procedure into a well-studied deformation flow, and we shall 
do it much in the spirit of a renormalization group flow ||14]| , [ 15] , [16| , [31| . As is 



common practice in the literature on the Ricci flow, we adopt the normalization that 
this flow is globally volume-preserving^. 

We prefer the Lagrangian point of view also for reasons of the following wider 
perspective. Suppose that the chosen hypersurface is a member of a given foliation of 
spacetime and, as an illustration, let the dynamical flow be geodesic. Then, the Einstein 
dynamics of, e.g., the spatial metric g a b is a first variation (in proper time t) of the metric 
in the direction of the extrinsic curvature tensor K a b, d t g a b = —2K a b- Evaluating the 
dynamical properties on spatial domains (geodesic balls) with fixed Lagrangian radius 
amounts to following the motion of the collection of fluid elements inside the ball along 
their spacetime trajectories, thus keeping the number of fluid elements (and hence their 
total mass) fixed during the evolution. This program has been carried out for the 
matter models 'irrotational dust', 'irrotational perfect fluid' and 'scalar field' in ||, J/J, 
fTT| . Complementary to this time-evolution picture (in the direction of the extrinsic 
curvature), we shall here investigate a d'Alembertian "virtual evolution" (first variation 
in the direction of the Ricci tensor TZ a b with variation parameter (3), where also in this 
case the material mass in the domain of averaging is kept constant. Our headline of 
forthcoming work is to setup a renormalized effective dynamics that includes spatial 
and temporal variation simultaneously. 



§ We have shown that a homothetic transformation of the variables would allow for a different 
normalization, e.g. such that the Ricci flow is globally mass-preserving. This, however, plays no 
significant role in our investigation of the regional averaging; the result is strictly equivalent. 
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1. Averaging and Scaling put into Perspective 

Addressing the issue of scaling in the Newtonian framework, one thinks of scaling 
properties of the dependent and independent variables under a transformation of the 
spatial Eulerian coordinates of the form x % \— > Lx' L (and a suitable scaling for the time- 
variable), where L is the parameter of the scale transformation (see, e.g., f43|). Such a 
transformation may be used as a dimensional analysis and for the purpose of constructing 
self-similar solutions of the Euler-Poisson system of equations. In general relativity we 
cannot have such a naive rescaling structure. Spacetime is a dynamic entity; as we shall 
see, rescaling of the form x l \— > Lx l must be replaced by a point-dependent functional 
rescaling, otherwise the scaled fields are strictly equivalent to the original fields. 

In order to characterize the correct conceptual framework for addressing averaging 
and scaling properties in relativistic cosmology, let us first remark that in general 
relativity we basically have one scaling variable related to the unit of distance: we 
can express the unit of time in terms of the unit of distance using the speed of light. 
Similarly, we can express the unit of mass through the unit of distance using Newton's 
gravitational constant. This remark implies that the scaling properties of Einstein's 
equations typically generate a mapping between distinct initial data sets, (in this sense 
we are basically dealing with a spatial renormalization group transformation). As 
we shall see below in detail, a rigorous characterization of an averaging procedure in 
relativistic cosmology is indeed strictly connected to the scaling geometry of the initial 
data set for Einstein's equations. 

Thus, it seems appropriate to recall some of the properties of the Arnowitt-Deser- 
Misner formulation of Einstein's equations. This essentially exploits the fact that a 
globally hyperbolic spacetime (M^ 4 \ g^) is diffeomorphic to R x E, where E is a 
three-dimensional manifold which we assume (for simplicity) to be closed. The explicit 
diffeomorphism cp : R x E — > is constructed by defining on a timelike vector 
field t the integral curves of which, $(t,p), define a map / : E — > according to 
p i — > f(p) = $(t,p), (for each given t G R). Consider, in such a setting, the dynamics 
of a self-gravitating distribution of matter. In the ADM formulation, the dynamics of 
such a distribution and the corresponding spacetime geometry (M^ ~ E x R,gr' 4 ') are 
described by the evolution of an initial data set: 

( E , g ab , K ab , q , J a ) (1) 

where the Riemannian metric gab and the triple of tensor fields (K ab , g, J a )^ are 
subjected to the Hamiltonian and divergence constraints: 

TZ + K 2 - K\K h a = IQtvGq + 2A , (2) 
V b K\ - V a K = 8nGJ a , (3) 

where K := g ab K ab is the trace of the extrinsic curvature K ab , 1Z := g ab 1Z ab is the trace 
of the intrinsic Ricci curvature of the metric g ab , and A is the cosmological constant. If 

|| Latin indices run through 1,2,3; we adopt the summation convention. The nabla operator denotes 
covariant derivative with respect to the 3-metric. The units are such that c = 1. 
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Figure 1. An admissible set of initial data is propagated along the time-like vector 
field t = niV+N, with the lapse function N scaling the unit normal n at the given point 
and the shift vector field N. The spacetime geometry is foliated into hypersurfaces S, 
(tfW - Ex JR,g (4) ), on which we study scaling properties of the initial data set. 



such a set of admissible data is propagated according to the evolutive part of Einstein's 
equations, then the symmetric tensor field K ab can be interpreted as the second 
fundamental form of the embedding of (£, g ab ) in the spacetime (M^ ~ £ x R, g^) 
resulting from the evolution of (E, g ab , K ab , g, J a ), whereas g and J a are, respectively, 
identified with the mass density and the momentum density of the material self- 
gravitating sources on (£,g a t,). Explicitly, the evolution equations associated with the 
initial data (£, g a b, K ab , g, J a ) are provided (in the absence of stresses, i.e., for a dust 
matter modeled) by 

?M = -2NK ab + £ M g ab , (4) 
-gf = - V a V b N + N [KK ab + TZ ab - (ArtGg + A)g ab ] + £ $K ab , (5) 

where N and N, respectively, denote the lapse function and the shift vector field 
associated with the mapping / : S — > (i.e., t = Nn + N, n being the future- 

pointing unit normal to the embedding £ ^h> (Fig. [I])). 

For the discussion of scaling properties and the averaging procedure associated 

We restrict our consideration to the simplest matter model, although almost all following 
investigations would apply to more general matter models. 
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with an admissible set of initial data (T,, g ab , K ab , g, J a ) for a cosmological spacetime 
(M^ 4 ) ~Sx R, gW) we start by characterizing explicitly a scale-dependent averaging 
for the empirical mass distribution g. 

1.1. Interlude: matter seen at different scales 

In order to characterize the scale over which we are smoothing the empirical mass 
distribution g, we need to study the distribution g by looking at its average behavior 
on regional domains (geodesic balls) in (E, g ab ) with different centers and radii. The 
idea is to move from the function g : E — > R + to an associated function, defined on 
E x R + , and which captures some aspect of the behavior of the given g on average, 
at different scales and locations. The simplest function of this type is provided by the 
regional volume average: 

1 



where p G E is a generic point, d[A g is the Riemannian volume element associated with 
(E, g a b), and B{p; r) denotes the geodesic ball at center p and radius r in (E, g ab ) (Fig. §), 
i.e., 

B(p;r):={qe(X,g ab ) : d g (p,q) < r } , (7) 
where d g (p,q) denotes the distance, in (E,^), between the point p and q. Note that if 

diam := sup {d g (p, q) :p,qE (E, g ab )} (8) 
denotes the diameter of (E,^), then as r — > diam, we get 

(q)b(p;t) — ► (*?>£ : = YT^y J ^ ( 9 ) 
at any point p G E. Conversely, if g : E — > R + is locally summable, then 

lin i (s) B (v.r) = Sip) , (10) 



for almost all points p G E. The passage from £> to (f?) B(jrr ) corresponds to replacing the 
position-dependent empirical distribution of matter in B(p;r) by a regionally uniform 
distribution (g) B ( p . r ) which is constant over the typical scale r. Note that for the total 
(material) mass contained in B(p;r), we get 

M{B{p;r)):= [ g dfx g = V (B(p;r)) (g) B{p;r) , (11) 

JB(n:r) 



and 



lim M(B(p;r)) = M(E) := g d/i g , (12) 

r—>diam J ^ 

where M(E) is the total (material) mass contained in (E,g a b). Our expectations in 
(@) B( P -r) are motivated by the fact that on Euclidean 3-space (R 3 , 5 ab ), if g is a bounded 
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**^B(p;r) 



Figure 2. The domain of averaging, centred on a given point p, is a (possibly 
non-convex) simply-connected geodesic ball domain contained in the hypersurface. 
The diagnostic parameter to explore the manifold is the radius r of the geodesic 
ball. It is mapped from the "master ball" with Euclidean geometry into the 
manifold by the exponential mapping. The shape of the ball's surface reflects the 
inhomogeneous geometry of the hypersurface at the scale r; the ball encloses a portion 
of the inhomogeneous matter distribution (marked by stars) interacting with the 
environmental distribution. 



function, then the regional average (g) B ( J) . r \ is a Lipschitz function on 
with the hyperbolic metric [41]: 

dx\ + dx\ + dx\ + dr 2 



endowed 



(13) 



In other words, 



(Q)B(p;r) ~ (0>B(«r) ^ C d h fa <?) 



(14) 



where C is a constant and dh(p, q) is the (hyperbolic) distance between p and q. Thus, 
the regional averages {g) B ^ p . r ^ do not oscillate too wildly as p and r vary, and the 
replacement of g by {(^) B ( p . r )} P eM 3 indeed provides an averaging of the original matter 
distribution over the length scale r. It is not obvious that such a nice behavior carries 
over to the Riemannian manifold (X, g a b). The point is that, even if the regional averages 
{(^)B( P r)}pes provide a controllable device of smoothing the matter distribution at the 
given scale r (see: Remark 1), they still depend in a sensible way on the geometry of 
the typical ball B(p; r) as we vary the averaging radius. In this connection we need to 
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Figure 3. Four ball domains mapped by the exponential mapping from the Lagrangian 
"master ball" into the Riemannian manifold are shown. Notice that the ball domains 
may overlap. The possible case in which they lie in topologically disconnected pieces 
of the manifold is excluded due to our choice of a maximal radius. 

understand how, as we rescale the domain B(p; r), the regional average {Q) B ( p . r ) depends 
on the underlying geometry of (£, g a b). The reasoning here is slightly delicate, so we go 
into a few details that require some geometric preliminaries. 

Geometrical evolution of geodesic ball domains 
Let us denote by 

exp p : T p £ — > £ 

(v,r) i — > exp p (rv) (15) 

the exponential mapping at p G (£, g^) , i.e., the map which to the vector rv G T p £ ~ M 3 
associates the point exp p (rv) G £ reached at "time" r G R + by the unique geodesic 
issued at p G £ with unit speed v G § 2 (1) . Let r > be such that exp p is defined on 
the Euclidean ball 

B E (0, r) := {y G T p £ ~ M 3 : < r} (16) 

and exp p : Be(0, r) — > B(p,r) C £ is a diffeomorphism onto its image. The largest 
radius r for which this is true, as p varies in £, is called the injectivity radius inj M of 
(£,<? o6 ) (Fig. |). 

Let B(p;r ) denote a given geodesic ball of radius r < inj M , and for vector fields 
X, Y, and Z, in (B(p;r), g ab ) let ^(X,F)Z = ^Vd^ c ^ d ^ fc 9 a be the corresponding 
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curvature tensor. Since r : (B(p;r), g ab ) — > R is a distance function (i.e., |Vr| = 1), the 



geometry of B(p; r) can be described by the following set of equations p8| : 

(V dr S)(X) + S 2 (X) = -TZ(X,d r )d r , (17) 
(C dr g)(X,Y)=2g(S(X),Y), (18) 
V 9r S = C 9r S, (19) 

where d r = Vr is the gradient of r, £g r g is the Lie derivative of the 3-metric g in 
the radial direction d r , the shape tensor[] S = V 2 r is the Hessian of r, and S 2 (X) 
abbreviates S 2 (X) = S a r .S V aXPd a . Such equations follow from the Gauss- Weingarten 
relations applied to study the r-constant slices U r := {exp p (rv) 6E:r = const.}, which 
are the images in (B(p;r), g ab ) of the standard Euclidean 2-spheres § 2 (r) C T p S ~ R 3 . 
The shape tensor S ab measures how the bidimensional metric g( 2 '(X,Y) induced on U r 
by the embedding in (B(p;r), g ab ) rescales as the radial distance r varies. If we denote 
by Pjj{W) :— W — g(W,d r )d r the tangential projection of a vector W G T P S onto the 
tangent space T p U to the surface U r at p, then together with (|Tj]), (0), ( JT9| ) we also 
get the tangential curvature equation: 

p u (n(x,Y)z) = 

G(gW(Y,Z)X-gW{X,Z)Y) - g^(S(Y), Z)S(X) - g^(S(X), Z)S{Y) , (20) 

where the vectors X and Y are tangent to U r , and where Q denotes the Gaussian 
curvature of (Ur^^). Further properties of the equations ( |T7| ) and ([LJJ) that we need 
are best seen by using polar geodesic coordinates. Recall that normal exponential 
coordinates at p are geometrically defined by 

exp; 1 : B(p; r) - T p £ ~ R 3 ; ex P ;\q) = {y*) , (21) 

where (y l ) are the Cartesian components of the velocity vector v G T P S characterizing 
the geodesic segment from p to q. Such coordinates are unique up to the chosen 
identification of T P S with IR 3 . Since we are dealing with a radial rescaling, for our 
purposes a suitable identification is the one associated with the use of polar coordinates 
in T p H ~ R 3 . We therefore introduce an orthonormal frame {ex, e-ii e z\ m T p T* such that 
ei := d r and with \e,2, e^\ an orthonormal frame on the unit 2-sphere §> 2 (1) C T p E. We 
can extend such vector fields radially to the whole T p £; we consider also the dual coframe 
9 2 , 8 3 associated with C2, e%. The introduction of such a polar coordinate system in T P S 
is independent of the metric g a b and thus is ideally suited for discussing geometrically 
the r-scaling properties of {Q) B u ) . r -\ (see, however, Remark 2). If we pull-back to T p S 
the metric g of B(p; r) C E, we get: 

g = dr 2 + g(e a , e p )9 a 6^, a, (3 = 2, 3, (22) 

where the components g a p := g(e a ,ep) = g^ 2 ' (r, cp) are functions of the polar 
coordinates (i?, ip) in T P E ~ R 3 , associated with the coframe 8 2 ,8 3 . Note that such 

+ The Hessian Sij = VyT" is the second fundamental form of the immersion U r =— > (B(p;r), g a b)- We 
use the equivalent characterization of shape tensor in order to avoid confusion with the standard second 
fundamental form of use in relativity. 
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a local representation of the metric holds throughout the local chart (B(p; r), exp" 1 ) 
and not just at p. In Cartesian coordinates in T p E ~ IR 3 one recovers the familiar 
expression: 

g ab = S ab - \n akbl {p)y k y l + O (M 3 ) . (23) 

In polar geodesic coordinates we have (x 1 := r, x 2 := x 3 := and g^l = g a p with 
2 < a, (3 < 3. Thus the equations (|17D , and (|TJ) take the explicit form: 

d r S a p = —S a rj S r 'i 3 — Tl a r p r ; (24) 
drgap = 2S\g v/3 , (25) 

where T^ a r p r denote the radial components of the curvature tensor. According to the 
tangential Gauss-Codazzi equation (ED) we can also write 

R 01 ^ — S (^s^ - fi^gpiu) - s^Su/j, + s a u S/3^ . (26) 

If we introduce the tangential components of the Ricci tensor according to 

7l a p := g kl n a m = K% r + - S^S + S a u S^ , (27) 

then we can rewrite the radial components of the curvature as 

Tl a r p r = Tl a p - Q5 a p + S a pS - S a u S u p , (28) 

and the first of equations (|25| ) becomes 

d r s% = -n a p + g5% - s a p s , (29) 

where 5* := S a gS a g is the rate of area expansion of g^ 2 \r, cp). Note that by taking the 
trace of both equations (|25| ) we get: 

d r S = -S a u S u a - Rtc(d r , d r ) , (30) 
g aP d r g afS = 2S , (31) 

where Ric(d r ,d r ) := 7Z l rir denotes the d r component of the Ricci tensor, (the first 
equation in (|30|) is nothing but the Jacobi operator coming from the second variation 
of the area associated with g( 2 )(r, <£>)). If we assume that the curvature 1Z(X, Y)Z 
is given, then, by fixing the (■#, y?) -dependence in the factorized metric (^2[), we can 
consider (|25|) as a system of decoupled ordinary differential equations describing the 
rescaling of the geometry of B(p; r) in terms of the one-parameter flow of immersions 
§ 2 (r) i — > (?7 r , (/ 2 )(r, y?)). Note in particular that the shape tensor matrix (S a g) is 
characterized, from the equation (p9[), as a functional of the Ricci curvature of the 
ambient manifold (E,g a b). 

Scaling properties of the metric 

In such a setting the equation d r g a p = 2S ri a g v /3 can be interpreted by saying that 
the metric rescales radially along the (curvature-dependent) shape tensor S%. In order 
to get a more explicit expression of such a radial rescaling, let us consider, at a fixed 
ro, a one-parameter family of surfaces (U r , g^ 2 \r, ip)) with parameter (r — tq) > 
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Figure 4. A foliation of the geodesic ball along the radial direction into a family of 
geodesic balls is shown. The metric on the boundary of the balls factorizes into radial 
and tangential components due to our choice of the Dreibein e\,e%, in the "master 
ball" . 

starting at the given surface (U ro , g^{r , (p)) and foliating the ball B(p;r) (Fig. [|). 
In a sufficiently small neighborhood of the initial surface (U ro , g^(r , </>)), we can 
write 

S%r) = S%r ) + [-K a p + - S%S] tq (r - r ) + 0((r - r ) 2 ) . (32) 
Inserting such an expression into d r g a p = IS^^^ we get 

^ r ~ r ^ ld r g a p(r) - d r g a/3 (r)\r ) = 

-2K aP {r Q ) + 2 [g5 aP - S a(3 S] ro + 0((r - r )) . (33) 

Note that the terms 2 [Q5 a p — S a pS} ro only depend on the geometry (intrinsic and 
extrinsic) of the surface (U ro , g^ 2 \r , </?)) and represent reaction tangential terms 
which work against the curvature of the ambient (E, gab). Thus, one may say that 
in a sufficiently small neighborhood of the initial surface (U ro , g( 2 \ro, ip)) the ambient 
geometry of (T,,g ab ) forces the metric to rescale radially in the direction of its Ricci 
tensor. This latter remark will turn out quite useful in understanding the geometric 
rationale behind the choice of a proper averaging procedure for the geometry of (£, g a b). 
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Scaling properties of the averaged density Geld 

Guided by such geometrical features of geodesic balls, let us go back to the study 
of the scaling properties of (g) B (p. r y To this end, for any r and s such that r + s < inj s , 
let us consider the one-parameter family of diffeomorphisms 

#.:(E,p)-(E,p) 

q = exp p [r,(d r , e 2 , e 3 )] i — ► H q (q) := exp p [(r q + s)(d r , e 2 , e 3 )] , (34) 

defined by flowing each point q G B(p; r) a distance s along the unique radial geodesic 
segment issued at p G E and passing through q. Let us remark that, for any r such that 
tq < r < inj A/ , we can formally write 

B(p; r) = H [r _ ro) B(p; r ). (35) 

Thus, in a neighborhood of B(p; r ) and for sufficiently small r, we get 



M(B(p;r)) :~- 



q du,„ 



B(p;r) 



B(p;r ) 



(36) 



where Hf_ r Jgdfj, g ) is the Riemannian measure obtained by pulling back (gdp, g ) under 
the action of if( r _ ro ). By differentiating ( |3~6"D with respect to r, we have: 



d 
dr 



M(B(p;r)) 



lim 

h->0 



[M(B(p; r + h))- M(B(p; r))] 



lim 



B(p;r ) (r-r )+h 



h 

[g dfig) - J B(p . ro) H* (r _ ro) (g dflg) 



h 



(37) 



Since Hi r _ ro \ +h = Hi r _ ro \ o H h , we can write the above expression as: 

f H (r-ro) [H hjgdUg) ~ (Q dfl g )] 

IB(p;r ) 



lim 



h 



lim 



U B(p;r) 



[H* h (gdfi g ) - (gd/x fl )] 



B(p;r) 



lim 



(38) 



from which it follows that (use Cd r dfi g = div(d r )dfx g and £a r £ = dg(d r ) = -§-g)'- 



-M(B(p;r)) = 
dr 

(C dr g + gdiv(d r )) d/i g 



B(p;r) 



B(p;r) 



B(p;r) 



9 1 ah 9 , , 



(39) 



where £g r and div(<9 r ) denote the Lie derivative along the vector field d r and its 
divergence, respectively, and where we have exploited the well-known expression for 
the Lie derivative of a volume density along d r : 











(40) 
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With these preliminary remarks along the way, it is straightforward to compute the 
total rate of variation with r of the regional average (p) B ( p . r ), since ( P5| ) implies: 



d , . Id 



dr (Q)B(p;r)-( dr Q 



B(p;r) 



+ K^L, - ^>^> (^L, ■ (4i) 

where (f) B (p. r \ denotes the volume average of / over the ball B(p;r). Explicitly, by 
exploiting (p5|), we get 

9 / 9 \ 

<0>B(p;r) = \fr e / B } + ( gS >B(p;r) ~ (&) B(p;r) ( S ) B(p;r) ■ ( 42 ) 

Thus, the regional average (^)_g(p- r ) f ee ^ s ^ ne fluctuations in the geometry as we vary the 
scale, fluctuations represented by the shape tensor terms 

(^>B(p;r) - (0)i?(p;r) ( S W;r) ( 43 ) 

governed by the curvature in B(p; r) according to (|30|) , and expressing a geometric non- 
commutativity between averaging over the ball B(p;r) and rescaling its size (Fig. 
Since the curvature varies both in the given B(p; r) and when we consider distinct base 
points p, the above remarks indicate that the regional averages (f?) B ( p . r ) are subjected to 
the accidents of the fluctuating geometry of (£, g ab ). In other words, there is no way of 
obtaining a proper smoothing of g without smoothing out at the same time the geometry 
of (£,5W>)- 

Non-commutativity of averaging and scaling 

It is straightforward to generalize the foregoing result for any smooth scalar function 
if) yielding a rule that summarizes the key result of this subsection. 

Proposition 1 (Commutation Rule) 

On a Riemannian hypersurface (S,^) volume-averaging on a geodesic ball B(p;r) and 
scaling ( directional derivative along the vector field -j^) of a scalar function ip are non- 
commuting operations, as can be expressed by the rule: 

s = 7f ; (s>B '-» = vmrT) ^ v{Bip ~- r » • (44) 

Note that this formula may also be read as follows: 
d d 

g^W)B(p;r) + (S)B(p;r)(^)B(p;r) = + Slf)) B (p;r) ■ (45) 



;(^) B{p;r) - (-K-^) B(p;r) = (Slf>) B(p;r) ~ (S) B(p;r){*P) B{p;r) ', 
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Figure 5. Volume-averaging scalar fields over the geodesic ball, and rescaling the 
radius of the ball are non-commuting operations. Note that a Lagrangian averaging 
as well as a rescaling to first order in the radius preserves the material mass content 
of the ball. 

1.2. Eulerian averaging and Lagrangian smoothing 

The use of the exponential mapping in discussing the geometry behind the regional 
averages {Q) B u,. r \ makes it clear that we are trying to measure how different the averages 
(^)b(pt) are f rom the standard average over Euclidean balls. In so doing we think of 
exp p : T P E — > E as maps from the fixed space Be{0,v) into the manifold (E,g a b). In 
this way we are implicitly trying to transfer information from the manifold (E,g a b) 
into domains of 1R 3 which we would like to be, as far as possible, independent of the 
accidental geometry of (E, g^) itself. Indeed, any averaging would be quite difficult to 
implement, if the reference model varies with the geometry to be averaged. This latter 
task is only partially accomplished by the exponential mapping, since the domain over 
which exp p : T p E — > E is a diffeomorphism depends on p and on the actual geometry 
of (E, #(,{,). A suitable alternative is to use harmonic coordinates in the ball B(p\r), a 
technique which is briefly discussed in Remark 2. 

We now go a step further by considering not just the given (E,g ab ), but rather a 
whole family of Riemannian manifolds. To start with let us remark that, if we fix the 
radius ro of the Euclidean ball Be(0; ro) C T P E and consider the family of exponential 
mappings, exp,^ : T p E — > (E, g a b(P)), associated with a corresponding one-parameter 
family of Riemannian metrics g a b(/3), < (3 < +oo, with g a b(P = 0) = g a b, then 
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B(p; r ) becomes a functional of the set of Riemannian structures associated with g a b{/3), 
< (3 < +oo, i.e., 

B(p; r ) i — > B p {p; r ) := exp (p >/3) [5 B (0; r )] . (46) 

In this way, instead of considering just a given geodesic ball B(p;r ), we can 
consider, as (3 varies, a family of geodesic balls Bp(p;r ), all with the same radius 
r but with distinct inner geometries g a b(P)- Since Bp =0 (p; r ) = B(p;r ), such balls 
can be thought of as being obtained from the given one B(p; r ) by a smooth continuous 
deformation of its original geometry. Under such deformation also (Q) B (p. ro ) becomes 
(3— dependent due to the functional dependence (g) B ( p . ro y 

The elementary but basic observation in order to take proper care of the geometrical 
fluctuations in (o) B ^ p . ro - ) is that the right member of ( |41~D has precisely the formal 
structure of the linearization ( i.e., of the variation) of the functional (g) B ( p . ro j in the 
direction of the deformed Riemannian metric -§g[g a b{(3)}, viz., 



d , , Id 



d(3 (e)B 0(p ;r o ) ( d p6 



+ \ (eg ab W)-^9a b W)) - \ (q) B()M (g ab (P)^9a b (P)) , (47) 

A \ °P I B (p;r o ) A \ °P I B fJ (p;r ) 

where the ball Be{0',Tq) is kept fixed while its image B(p; r ) is deformed according to 
the flow of metrics g a b(P), < (3 < oo (Fig. ||). 

In a rather obvious sense, (f47D represents the active interpretation corresponding 
to the Eulerian passive view associated with the ball variation B(p;r Q ) — > B(p;r). In 
other words, we are here dealing with the (geometrical) Lagrangian point of view of 
following a fluid domain in its deformation, where the fluid particles here are the points 
of B(p;ro) suitably labelled. This latter remark suggests that in order to optimize the 
averaging procedure associated with the regional average (g) B ^ p . r y instead of studying 
its scaling behavior as r increases, and consequently be subjected to the accidents 
of the fluctuating geometry of (B(p;r), g a b), we may keep fixed the domain Be(0; r ) 
(setting the scale over which we are averaging) and rescale the geometry inside its image 
B{p; r ) under the exponential map, according to a suitable flow of metrics g a b —> g a b{(3), 
< (3 < oo. Correspondingly, also the average matter density will be forced to rescale 
(Q)B(p;r ) — ► (s)Bg(p;r ){l3), and if we are able to choose the flow g a b — > g a b(/3) in such 
a way that the local inhomogeneities of the original geometry of (E, <7 a &) are smoothly 
eliminated, then the regional averages (g) B come closer and closer to represent a 
matter averaging over a homogeneous geometry. 

According to the r— scaling properties of the metric described by (0), a natural 
candidate for such a Lagrangian flow is the deformation generated by the Ricci tensor 
of the metric, a deformation flow that is strongly reminiscent of the Ricci flow on a 
Riemannian manifold (S, <7 a b)p7|, fl~3H , [f28fl : 
d 

g^9ab(v) = -2^(77) ; g ab (r] = 0) = g ab , (48) 
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Figure 6. Averaging over geodesic ball domains, while scaling the geodesic radius 
is conceptually (and formally) equivalent with their Lagrangian deformation for fixed 
geodesic radius. The transformation of a Lagrangian ball is illustrated with respect to 
r— scaling and (3— deformation. 
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Figure 7. When it is global, the Ricci flow acts in the space of all possible Riemannian 
geometries by smoothing the geometry in a controllable fashion: here we depict the case 
of a three-sphere whose maxima of the curvature inhomogeneities are monotonically 
decreasing. If the original three-sphere has positive Ricci curvature (i.e., is not wildly 
inhomogeneous), then it undergoes a metamorphosis into the standard round three- 
sphere. Globally, if large inhomogeneities arc initially present, the Ricci flow may 
experience singularities, the flow solution may feature a bifurcation and consecutive 
regions may be 'pinched off' in the limit where they are smoothed-out. This implies a 
topology change of regional domains. 
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studied by Richard Hamilton and his co-workers in connection with an analytic attempt 
to proving Thurston's geometrization conjecture. As is well-known, the flow ([48|) 
is weakly-parabolic, and it is always solvable for sufficiently small rj. Obviously, it 
preserves any symmetries of g a bij] = 0) = g a b- (The Ricci flow preserves the isometry 
group of (E, g a b)-) It can be viewed as a "heat equation" diffusing Riemannian curvature 
(Fig- 0). 

Global normalization of the Ricci flow 



The flow (f48D may be reparametrized rj — ► (3 by an ^-dependent rescaling and by 
an ^-dependent homothety g a b{v) ~ * 9ab{(3) so as to preserve the total volume of the 
manifold (E, <7 &(/3)) ='■ E^. F° r this end one has to introduce a suitably normalized 
flow: 

^g ah {P) = -2K ah {(3)+A{P)g ah {(3) ; g ah {P = 0) = g ab ; (49) 

the normalization factor A(/3) is determined in such a way that the volume of E^ 
associated with the metric g ab (P) does not change under deformation, <9/<9/3Ve„(/3) = 0. 
In place of P8| ) we then get the standard volume-preserving Ricci flow that is usually 
studied in the mathematical literature P7| , 



—g ab {f3) = -2n a b{(3) + -(n{(3)) g ab {(3) ; g ab (P = 0) = 9ab , (50) 

OA 3 \ / T, p 

whose global solutions (if attained) are constant curvature metrics: 

9ab ■= lim 9ab(/3) ■ (51) 

/3-»+oo 

Since the normalization factor is spatially constant, it will not enter as a fluctuating 
quantity in our equations for the regional averages. If we would, e.g., normalize the Ricci 
flow such that it preserves the global mass, this would not change the statements on 
regional averages, where we shall require the preservation of the regional mass. We 
emphasize that such a normalization is a technical choice, mathematically needed in 
order to be able to compare the distinct regional averagings carried out with respect to 
balls with different centers. Note that, even if we only wish to smooth the hypersurface 
E^ on regions of (Euclidean) radius r , (i.e., B E (0;r )), their representatives Bp(p;r ) 
are to be considered for distinct centers, say pj] in other words, we can average over 
Bp(pi] r ), Bp(p k ; r ), where {pi, ...,Pk} is a set of points suitably scattered over the 
manifold E^ (compare Figure 3). As a matter of fact, all our final results factor out the 
global volume average and refer only to the average with respect to the regional ball. It 
must also be stressed that, to our knowledge, there is not yet a mathematically correct 
way for implementing a Ricci flow that just works for an open region (such as a ball) of a 
3— manifold (one needs suitable boundary conditions, on the spherical boundary of the 
ball, controlling the flow of curvature in and/or out of the ball). Since the total volume 
is not a physically observable quantity, we take care of working out results which are 
independent of the volume constraint. 
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2. Averaging and Scaling put into Practice 



2.1. Smoothing the metric 

Let us now come to the strategy for the optimal choice of the smoothing flow for the 
metric g a b — ► 9ab(P)- As outlined in the previous section, when dealing with regional 
averages, a suitably normalized Ricci flow taking care of the metric comes naturally to 
the fore. 

To set notation, let us again write the volume-preserving Ricci flow equations in the 
form: 

^g a k(j3) = -2fcab(P) + lgab(P)(K(J3))x f > I 9ab((3 = 0) = g ab , (52) 

where < (3 < +oo is the deformation parameter. As already stressed, on any compact 
Riemannian manifold (E, g a b) we have Hamilton's theorem according to which a local 
solution to such a flow with g a b(/3 = 0) = g a b always exists for (3 sufficiently small E7 



18] ; for an introduction to such problematics, with many self-explanatory diagrams, 
14f . The study of the existence and properties of global solutions g a b — > g a b(P), 



sec 



< (3 < oo, is much more difficult to establish, and is an active field of research (see 
the recent review by R. Hamilton [25], and [|12|]. In particular, if the initial metric 



(Y,,g a b) has positive Ricci curvature, then the solution g a b — > g a b{(3) to (|52j) exists for 
all (3, and it converges exponentially fast, as (3 — > oo to a constant positive sectional 
curvature metric (E, </ o6 ), (forcing E to be a space form diffeomorphic to the 3— sphere 
§ 3 , possibly quotiented by a finite group F of isometries). Other examples of a global 
Ricci flow are provided by those flows that evolve from locally homogeneous metrics 
(E, <7 Q &). In particular, the eight distinct homogeneous geometries existing in dimension 



n = 3 have been analyzed in detail resulting in a non-singular global flow | 32fl ; for an 



example originating from relativity see |]15 



Since, according to Thurston's geometrization conjecture (Fig. ||), every closed 3- 
manifold can be decomposed into pieces admitting one of the eight geometric structures 
mentioned above, it is clear that the global Ricci flow may play a distingushed role 
in such a conjecture. While such a role clearly motivates the mathematical interest in 
([52]), it also provides a strong argument in favour of (0) as the natural smoothing flow 
in a regional cosmological averaging procedure. As a matter of fact, the possibility of 
decomposing a 3-manifold into pieces endowed with a locally homogeneous geometry is 
particularly appealing to relativistic cosmology, where any such piece may be thought 
of as representing the regional average of a sufficiently homogeneous portion of the 
Universe. In such a framework it is suggestive to put the following result by Richard 
Hamilton into perspective p9l (see also |p~2|): 



Regional Averaging and Scaling in Relativistic Cosmology 



22 




Figure 8. Thurston's conjecture as viewed by Richard Hamilton is illustrated. The 
global Ricci flow may have singularities resulting in an effective disconnection of the 
manifold into "nice" pieces. Thurston's conjecture assumes that the original manifold 
can be cut into pieces which themselves have elementary Ricmannian structures 
(associated with the eight "Thurston geometries"). The figure shows three such pieces 
in a two-dimensional rendering which would correspond to pieces of the three possible 
FLRW space sections in three dimensions: the 3— sphere, a piece of a 3— cylinder, or a 
piece of a 3— hyperboloid. 
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Theorem (Hamilton) 

If the closed 3-manifold (E,g ab ) admits a non-singular solution g ab — > g a b(P) to (p^) 
for all < (3 < oo, with uniformly bounded sectional curvature, then (E, ~g ab ) can be 
decomposed into pieces admitting one of the following locally homogeneous geometries: 



U 

(ii 
(iii 

(iv 

fv 



(E, g ab ) is a Seifert fibered space. 

(E, g ab ) is a spherical space form § 3 /T. 

(E, ~g ab ) is a flat manifold. 

(E, ~g ab ) is a constant negative sectional curvature manifold. 

(E, g ab ) is the union (along incompressible tori) of finite-volume constant negative 
sectional curvature manifolds and Seifert fibered spaces. 



The locally homogeneous geometries («)•■• (v) (in particular (ii), (iii), and (iv)) are 
exactly the geometries after which, we believe, the Universe can be regionally modelled 
after all accidental inhomogeneities are ideally ironed out. Thus, Hamilton's theorem 
strongly advocates the basic role of the Ricci flow deformation as a natural mean for 
averaging locally inhomogeneous 3-geometries in a cosmological setting. 

On the cosmological stage, however, we need also to discuss how such a geometrical 
smoothing flow interacts with the actual distribution of matter. 

Combining the Ricci How and the material mass how 

Our basic idea behind the regional averages (f?) B ( p . r ) is that they replace the local 
accidental distribution of matter described by q. Also, instead of considering just a 
given geodesic ball B(p;r ), we are considering, as (3 varies, a family of geodesic balls 
Bp(p;ro), all with the same radius r but with distinct inner geometries g a b(P)- Note 
that r plays here the role of a distance cut-off: it is the typical scale over which we 
want to smooth the empirical mass distribution. We have to keep track of such a scale 
in our setup. 

The key idea for obtaining a deformation equation for the regionally averaged mass 
density field is to fix, besides the scale ro of the Lagrangian ball, also the mass content of 
the ball under consideration during the deformation process (Fig. ^). It is a physically 
sensible idea to concentrate on mass scales, since they are only directly related to so- 
called "comoving" length scales or volumes, respectively [[37], if the Eulerian ball is also 



Euclidean and the inhomogeneites are set up relative to a global reference flow providing 
the global comoving coordinates. 

In order to obtain a deformation equation for the regionally averaged mass density 
field in the chosen regions Bp(pj \ ro), we require 

d , , , „ , xx d 
_M( B ,fe;r ))=:-. 



— a M (B s {p,; r„)) =: — M Bf = (53) 
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Figure 9. The conservation of the total material mass contained within geodesic ball 
domains guarantees the comparison of average properties on cosmological mass-scales. 

along the solution of (j52|). Since (the abbreviation Bp := Bp{pj]ro) is used hereafter) 
<ff>B, 03) = Ti~Ta\ . (54) 



where Vb denotes the (for each (3 different) volume of the ball in question, it is 
straightforward to check that we can equivalently rewrite (Ek3) as: 

8 In (e)B > W I a i n ^^» Q (55) 

We realize the local rate of volume deformation through the Ricci deformation flow fl52p 
by exploiting the relation: 

Tpy^^^ = l^W^gM = W)h, - n{p) , (56) 

where 7?.(/3) := g ab (f3)1Z a b((3) is the scalar curvature associated with g a b(P)- From this 
one easily computes that 

J_ ln v^O = _ {nm ^ + W)) ^ (57) 



d(3 V t 



Thus we get that the condition of the conservation of the matter content of Bp(pj;ro) 
(|53|) provides the following /3-evolution for the regionally averaged mass density: 



d_ {q)b b iP) 



= W)) B -{K{(3))„ . (58) 
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So far we have established an analogy between matter averaging on different scales 
and geometrical deformation induced by a suitable flow of metrics. With regard to the 
constraint equations of general relativity we have to guarantee that this flow of metrics 
is compatible with the constraints. Before we turn to the problem of smoothing the 
second fundamental form, it is necessary but also illuminating in what follows to study 
stability properties of the Ricci flow. 



2.2. Stability of the Ricci flow 

Associated with the Ricci flow we need to discuss also the properties of the corresponding 
linearized flow. Roughly speaking, such a necessity comes about since we may be asked 
what happens to our flow, if the original metric g a b, deformed according to g a b — > g a b{(3), 
is slightly perturbed (which is a good question, since the flow is actually perturbed in 
a dynamical situation, e.g., in the direction of the extrinsic curvature tensor in time 
e = t): 

9ab -» 9ab = 9ab + ^Kb , (59) 

where h ab is a symmetric bilinear form, and e a small parameter. For notational 
convenience, let us set 

2 
3* 

so that the Ricci flow can be compactly written as 



RF[9ab(P)} ■= -2K ab (P) + -g ab {(3) W))^ , 



(60) 



-^gab(P) = RFM/3)] 



9ab(P = 0) = g a b ■ 



(61) 



It is easily checked that, if we perturb the initial metric g a b(P = 0) = g a b, and evolve 
the perturbed metric 'gab according to (p2|) , then we get 

g ab -> gobiP) = 9ab{(3) + ehab{[3) + 0{e 2 ) , (62) 

where h a b — > h a b(P) is a linear flow solution of 
d. _ d 



o h ab ((3) := -RF[g ab ((3) + eh ab ((3)) 



de 



e=0 



0phab((3) = DRF[g ab [ 

h ab (P = 0) = h ab ■ (63) 
Explicitly, we obtain (dropping the explicit /3-dependence for notational ease), 



d 2 



+ -^9ab 



-Arh 



[Kh ik g ik )^ 



L lb ab 



5* 5 h 



2 ( n K 

l;( tTh )9 



'h ik g ik )„ - {K ah h ab \ 



ab 



where 



-V s V s h ab + HasK + TZ bs h s a - 2TZ asb th st 



(64) 



(65) 
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is the Lichnerowicz-deRham Laplacian on bilinear forms, the operator 5 is (minus) 
the divergence, and 5* is its formal L 2 -adjoint ( i.e., 1/2 the Lie derivative operator). 
All such operators are considered with respect to the (3- varying metric g ab (P) of the 
unperturbed Ricci flow (0). It is clear from its explicit expression that flBll) takes quite 
a simpler form, if we restrict our attention to traceless perturbations h ab , i.e., if 

h ab g ab = . (66) 

It is then verified that such a condition, if it holds initially (for f3 — 0), it holds 
for each value of /3 for which the flows (|5^) and (|64l) are defined, and we get the much 
simpler result: 
d 2 2 

■gp h ab = g (ft) E(9 Kb - <"^ a ^a&) S/3 - A L /i a6 + £„ V k h . k g ab . (67) 



For a given Ricci flow g ab — > g ab (f3), Eq. ( p7|) defines a linear (weakly) parabolic initial 
value problem (the strict parabolicity is broken by the Lie derivative termQ: 

£-vH Hk 9ab = -V a (V k h kb ) - W b iy k h ka ) , (68) 

associated to the infinitesimal equivariance under the diffeomorphisms group Diff(S)). 
Given the initial (traceless) perturbation, h ab (j3 = 0) = h ab , the solution of (|67D always 
exists and is unique, and it represents an infinitesimal deformation of metrics connecting 
two neighboring flows of metrics g ab — > g ab (P) and g ab — > g a b(P)- Since h ab ((3)g ab (f3) = 0, 
both flows have the same j3 -dependent volume element, i.e., 

vW) = VW) , (69) 

and thus the same /3-dependent average density (g) B . It is also important to remark 
that the solution of ( |B"TD corresponding to the trivial initial datum (a conformal Lie 
derivative term £^g ab ), 

2 

h ab ((3 = 0) = £^g ab ■= V a w b + V b u; a - -g a bV c w c , (70) 

where w b is a smooth (/^-independent) vector field on S, is provided by 

h ab {f3) = £^((3) := V{(3) a w b + V((3) b w a - ^g ab (f3)V((3) c w c , (71) 

where the /5-dependence is only through the flow of metrics g ab — > g ab (/3) and the 
associated connection V a (/?). Such a property is simply a consequence of the Diff(E) 
equivariance of the Ricci flow (Fig. [ToD . By exploiting such a result, it is possible to prove 
an important factorization theorem ^5) for the structure of the solution h ab — > h ab (f3) 
of (|67D, which will prove invaluable throughout the rest of the paper, viz., 



* Note that £ abbreviates the conformal (i.e., trace-free) Lie derivative, which should not be confused 
with the Lie derivative denoted by C. 
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Figure 10. The stability analysis of the Ricci flow under symmetric perturbations 
allows us to determine the smoothing flow for the second fundamental form. It can be 
viewed to represent an infinitesimal deformation of metrics connecting two neighboring 
flows of metrics. The flow of this deformation itself can be shown to be (weakly) 
parabolic (like the Ricci flow for the metric) up to the infinitesimal equivariance of the 
Ricci flow under the diffcomorphisms group Diff(S) ( "Diff- Warping" ) . 
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Proposition 2 (Lott) 

If hab habifi) is the flow solution of corresponding to the initial (traceless) datum 
h a b{ft = 0) = hat,, then it can always be factorized according to 

Kb(ft) = K b (ft) + £ m9a M , (72) 

where the bilinear form is the solution of the parabolic initial value problem 

^Kb = -Al/4 + 1 (7^ Kb - fa* {n ab K b ) nfi) , (73) 

with h* ab (ft = 0) = h ab , and where the (now ft— dependent) vector field v a (ft) is the flow 
solution of 

^v a (ft) = -V C K C ((3) , v a (ft = 0) = . (74) 

It is appropriate at this point to recall a few relevant facts concerning the geometry 
behind the structure of the solutions of fl67|) . It follows from the above proposition 
that, as ft — > oo, h a b(ft) may either approach a (conformal) Lie derivative term £^g a b, 
or a non- vanishing deformation tensor h* ab . This latter non-trivial deformation is only 
present, if the corresponding Ricci flow g a b — » 9 ab(ft) approaches an Einstein metric 
on E which is not isolated (for instance flat tori). In such a case, there is a finite- 
dimensional manifold of such Einstein metrics, and the non-trivial h* ab simply represents 
infinitesimal deformations connecting two infinitesimally neighboring Einstein metrics 
on E. As is known, the round metric ~g ab on the three-sphere § 3 is isolated in the sense 
that there are not volume-preserving infinitesimal deformations of ~g ab mapping it to 
another inequivalent constant curvature metric ~g' ab . In this latter case, {i.e., for isolated 
constant curvature metrics), as ft — > oo, h a t,(ft) must necessarily approach a (conformal) 
Lie derivative term £$~g ab - 

2.3. Smoothing the second fundamental form 

The properties of the linearized Ricci flow for a traceless metric perturbation h ab (ft) 
naturally put to the fore an explicit way for averaging the part of initial data 
(E, g a b, K a b, £?, J a ) related to the second fundamental form K a t,. One may contend 
that since K a b carries information on the way (£, <jr a &) is embedded in the spacetime 
(M^ 4 ) ~ E x R, g^), one should devise some way of deforming K ab which is independent 
of the flow of metrics, since this latter flow only depends on the intrinsic geometry of 
(E, gab). However, the very geometrical meaning of K ab shows that such a point of view 
is not correct. According to the evolutive part of the Einstein equations we have: 

Kab = ~^N^r + ^N c ^ 9ab - (75) 

Thus, we can write 

gabif) =g ab -t[2NK ab - C$g ab ] + 0(t 2 ) , (76) 

which clearly shows that 2NK ab has the natural meaning of the deformation tensor 
connecting two neighboring Riemannian metrics. If the Ricci flow is the chosen averaging 
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procedure for deforming the metric in the initial data set (E, g ab , K ab , g, J a ), then the 
stability (under small perturbations) of such a deformation procedure requires that 
NK ab must necessarily be deformed according to the linearized flow (^ ); there are no 
other consistent choices. 

The only freedom we have concerns which (algebraically independent) part of K ab 
we want to deform according to (0). From the properties of this latter flow it follows 
that a smart choice would be to leave undeformed the trace part of K ab , and deform 
only its traceless part (i.e., the associated shear tensor a a b '■= K ab — \g a bK). 

Actually, the geometrical meaning of K ab epitomized by (|76|) suggests as a more 
natural choice that we deform the trace-free distortion tensor (Fig. |Tl|): 

V ab := 2Na ab - £ $g ab 

= 2Na ab - \V a N b + V b N a - ^g ab V c N^j . (77) 

In order to /3-deform T> ab in a way consistent with the Ricci flow, let us start observing 
that, with respect to the metric g ab , we can always decompose the given V ab according 
to 

V ab = V ±ab + V mb (78) 
where V± ab is the divergence-free part, 

V a V± ab = , (79) 
and V mb is the longitudinal part, 

T^mb = £w9ab , (80) 

generated by the vector field w as a solution of the (elliptic) partial differential equation: 

V a (£n!9ab) = V a \V a W b + V b W a - = \/ a V ab . (81) 



Observe further that, according to (|76|), T>± ab represents the part of the distortion 
tensor that deforms g ab into a nearby distinct Riemannian structure [g a b(t)} = g ab — 
D± ab t + 0(t 2 ), whereas T> uab simply generates an infinitesimal Diff(S) reparametrization 

9ab —> 9ab ~ t [£w9ab ~ £>]y9ab} + 0(t 2 ). 

With all this in mind, we can apply Lott's factorization theorem and generate a 
natural smoothing flow for V ab , [3 — ► V ab (f3) as a solution of 

§pV ab (P) = | (K)^ V ab (P) - 2 -g ab (TZ ab V ab ((3))^ 

- A L V ab ((3) + £„ vkVim g ab ■ (82) 

by setting V ab (f3) = P* 6 (/3)+ £vg a b(P), where T>* b (/3) evolves according to the parabolic 
partial differential equation 

^fl^) = + I (K)^ V; b (P) - \g ab (n ab V* ab (P)) sw , (83) 
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Figure 11. Smoothing the second fundamental form we have the freedom to choose 
which of the algebraically independent parts of the second fundamental form we wish 
to smooth out. The trace-free distortion tensor is the natural choice in view of the 
evolutive part of Einstein's equations. In the embedding shown in the figure we have 
assigned the coordinate time to the perturbation parameter e, and the perturbing 
symmetric bilinear form was associated with the second fundamental form. 
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with T>* b (j3 = 0) = V ab , and where the (3— dependent vector field v a {(3) is the flow 
solution of 

= -V c V* ac ((3) , v a ((3 = 0) = . (84) 

Corresponding to a global solution (3 — > g a b(P) of the Ricci flow, that evolves towards 
an isolated constant curvature metric g ab = lim^ +00 g ab (j3), the (unique) solution 
P -> (V*ab(/3),v a (P)) o f (HD, dH) is such that V ±ab = \im^ +00 V ±ab ((3) = 0, (V ±ab 
will be different from zero only if the constant curvature metric ~g ab is not isolated), and 
V ab reduces to a pure longitudinal shear, i.e., 

V ab = lim V ab {(3) = £0 ab , (85) 

where 

v a = lim v a ((3) . (86) 

0— *+oo 

Note that since we can always add to T> ab ((3) the trivial solution £^g a b(P) of (0)> (where 
the vector field w does not depend on /5), by suitably choosing the (/3-independent) 
vector field w, we can always assume that the conformal Lie derivative term £ ^g ab 
provided by flSB]) comprises all the longitudinal shear present in T> ab . 



2.4- The choice of a smoothing reference frame 

It is important to stress that the averaging of the trace-free distortion tensor described 
above implies that we have made a consistent selection for the lapse function N, the 
shift vector field iVf,, and the rate of volume expansion K associated with the given 
initial data set (S, g ab , K ab , g, J a ). As is well-known, such (N,N b ,K) are kinematical 
quantities pertaining to the choice of the initial slice (£, g ab , V ab ). They are related also 
to the choice of the foliation H t in a suitably small (time) neighborhood of £ in the 
spacetime resulting from the time-evolution of the data (E, g ab , K ab , g, J a ). Explicity, 
we have 

d 

—K = -A g N + N[R + K 2 - 3(4nGp + A)] + N a V a K . (87) 

This is basically Raychaudhuri's equation. The rationale underlying the choice of a 
proper reference frame for carrying out a sensible (regional) averaging, is that (^) should 
reflect the achieved regularity in the averaged geometry and matter fields, without 
introducing frame-dependent inhomogeneities and anisotropies. This implies that, as 
we /5-deform the local inhomogeneities and anisotropies of (£, g ab , V ab ), we need also 
to eliminate the inhomogeneous artifacts due to the choice of the slicing associated 
with (E, N, N b , K). To give an example, it is often argued that a good candidate 
for supporting an "almost-Friedmannian" initial data set (E, g ab , K ab ) is the surface of 
constant matter density in the cosmic fluid frame, or the surface of constant 4-velocity 
potential (0, JEJ, (§, [§], 0, [0, 0). The point is that in such a slicing a slightly 



perturbed model features an almost constant lapse function, since the instantaneous 
acceleration for such a frame of reference is observationally quite small. The observed 
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expansion and shear are also small. In line of principle these remarks suggest that a set 
(E, N, Nb, K) characterized in this way is, in the observed domain, the most suitable 
one for implementing an averaging procedure. At any rate, even for such a natural and 
almost optimal choice of (E, N, N b , K), we still have the issue of how to consistently get 
rid of the (residual) frame fluctuations present in N, Nb, and K. In order to eliminate 
such fluctuations, both the lapse function N(0), the shift vector field N a {(3), and the rate 
of volume expansion K((3) are to be considered as explicitly (3— dependent, (taking the 
given values (N, Nb, K) for (3 = 0). In other words, we need to consider on (E, g a b, V a b) 
a one-parameter (/?) family of (instantaneous) frames of reference (E, N(f3), Nb(/3)) and 
devise a way for characterizing their /3-evolution. The rationale is to end up in a 
frame which, for (3 — > oo, reflects, as much as possible, the homogeneous and isotropic 
properties of the geometry resulting from the (3 -evolution of (E, g a b, V a b). Such a frame 
will correspond to the standard Friedmannian scenario of use in cosmology. 

A word of caution is mandatory here. In mathematical cosmology it is often the 
case that the choice of foliation is strictly connected to the structure of the constraint 
equations via the use of K as the variable conjugated to time (York's extrinsic time). 
Spatially constant K is then a rather popular choice. However, the rate of volume 
expansion plays a distinguished dynamical role in Friedmannian cosmology, and for our 
purposes it would be quite detrimental to use K as the variable selecting the hypersurface 
E carrying the data to be averaged. In line of principle, the structure of the Ricci flow 
suggests that one should pick up a E such that (E, g ab ) admits a global < (3 < oo Ricci 
flow, for instance this is the case if the Ricci tensor of (E,g a b) is positive or (hopefully!) 
not too wildly oscillating. However, this is again something which is under control 
of mathematics, but not acceptable for the physical situation; it cannot be used as a 
viable selection criterium. The best we can hope for is to assume that E is ideally chosen 
among the t-constant slices of a global frame of reference where, for scales sufficiently 
larger than the relevant averaging scale, the distant galaxies appear to recede radially. 
And, as argued before, such a E may be appropriately realized, e.g., by the surface of 
constant matter density in the cosmic fluid frame. 

The lapse function 

We postpone to the next section the issue of the /3-evolution of the rate of volume 
expansion K, since this latter is strictly connected to the regional averaging of the 
constraints. On the other hand, for the (3 -evolution of the kinematical variables 
(N, N b ) there are physically sound choices which are directly suggested by the nature 



of the Raychaudhuri equation, and from the structure of the equations fl83|), and (|84|). 
To start with, let us observe that since the leading lapse-dependent inhomogeneity 
inducing term in (^7|) is A g N it is natural to smooth the lapse function iV by diffusing 
its inhomogeneities by means of the scalar heat equation: 

-^N(f3) = A m N(f3) ; N{fi = 0) = N , 
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where A g /g) is the Laplacian of the metric g a b((3) evolved by the Ricci flow. This is 
basically the harmonic map flow for the map N((3) : E — ► R. Note also that this 
is a non-uniformly parabolic initial value problem, because A g (g) depends on the (3- 
varying metric g a b((3)- Let us assume that, for (3 = 0, the given lapse function is such 
that < e < N < Ci, \VN((3 = 0)| < C 2 , where C\ and C 2 are suitable constants. 
In other words, the (instantaneous) acceleration V a In N of the frame associated with 
the chosen slice E is assumed to be bounded. Then, according to the maximum 
principle for © (e.g. see 0), we have < e < N((3) < C u \VN(/3)\ < C 2 for 
(3 > 0. If the Ricci flow exists for all (3 > 0, then by looking at the /3-evolution of 
\VN((3)\ 2 = g ab ((3)V a N ((3)V bN ((3) , a long but elementary computation provides: 

||V7V (/3 )| 2 

= A sOT (|ViV(/3)| 2 )-2|V 2 iV(/3)| 2 + V°iVV i ']V^9«i(/?) + 2K al ,J , (89) 

which shows that, if we deform the metric along the Ricci flow, then we have: 

A \VN((3)\ 2 = A m (\VN(P)\ 2 ) - 2 |V 2 iV(/3)| 2 - | |ViV(/3)| 2 (fc) E/j . (90) 

It follows that the maximum of \VN((3)\ 2 is weakly monotonically decreasing as /3 —> oo; 
(the apparent condition for this to occur is (TZ) S > 0, but this condition is actually not 
necessary, as can be seen by rescaling the (3— variable to the standard unnormalized Ricci 
flow). Parabolic theory shows that, if the Ricci flow is global (with uniformly bounded 
curvature), then also (3 — > N(/3) exists for all (3 and \VN((3) \ — ► as (3 — > oo. Physically, 
such an averaging procedure implies that the frame instantaneous acceleration associated 
with the chosen slice S is smoothly damped, and up to a normalization we can assume 
that as (3 — * oo, N(f3) — > 1 uniformly. Note that if one identifies E with a surface 
of constant matter density or constant 4— velocity potential, respectively, in the frame 
of reference comoving with the flow lines of the (irrotational) cosmic fluid, then an 
inhomogeneous N can be directly related to the instantaneous acceleration of the fluid 
particles on E, thus N((3) — ► 1 as (3 — > oo implies that the fluid particles are, on E, in 
free fall. 



The shift vector and the matter current density 

Since N a generates local diffeomorphisms in the hypersurface E, a procedure for 
smoothly averaging out the shift vector field iV a cannot be disentangled from the 
averaging of the matter current density J a and from the fact that the averaging of 
the distortion tensor V a b(f3) generates, according to Lott's factorization theorem, a 
geometrically induced shift v a ((3) (see (0)). Along the same lines of the above remarks 
one may tentatively assume that iV a can be related to the local 3-velocities of the fluid 
particles on E. On a sufficiently large cosmic domain we may put ourselves initially on 
E (i.e., for (3 = 0) into a frame with vanishing shift: N a = = J a . Naively, one would 
expect that such a situation will persist as we deform the data (E, g a b, V a b). However, 
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as already stressed, there is a diffeomorphism warp generated by the linearized Ricci 
flow, which, if not properly taken into account, will manifest itself as a current J a of the 
matter on the smoothed data (E, g a b(ft), T^ab(ft)) when ft — > oo. The physical origin of 
such diff-induced matter current is rather easy to understand. Roughly speaking, what 
happens is that the points of the manifold E are moved around as curvature bumps 
are ironed out, (the generator of such a motion is the gradient in scalar curvature 
V a TZ; these are basically the Diff-solitons familiar in the Ricci flow theory). From a 
Lagrangian point of view, this motion is transferred to the fluid particles labelled by 
the corresponding points of E. 

Thus, in order to consistently compensate for such a Ricci flow induced warping, we 
must assume that for ft > we have a non- vanishing shift N a (ft) and a corresponding 
matter current density J a (ft). The natural choice that suggests itself is to introduce a 
ft -dependent shift ft — > N a (ft), according to 

N a (J3f)=v a -v a (J3) , (91) 

where v a (ft) is the flow solution of (|84l) , and v a is its ft — ► oo limiting value given by 



(36|). Note that, according to such a choice, 

N a (ft = 0)=v a , N a {ft oo) = , (92) 

so that the shift N a (ft) exactly balances the longitudinal shear £ v g a b(ft) generated by the 
(linearized) Ricci flow for < ft < oo. In other words, there is an optimal choice v a for 
the 3-velocity field of the instantaneous observers on E which allow them to isotropically 
smooth the data (E, g a b, T> a b)- If we, e.g., identify E with a surface of constant average 
matter density (g) B , then the Ricci flow induced shift ft — > N a (ft) gives rise to a fictious 
material convection current density ft — > J a (ft) given by 

J a ((3):=(g)^N a (ft) , (93) 

where (g) B is the averaged matter density over the region of interest, (see (|58D ). Observe 
that J a (ft) — > uniformly as ft — > oo. 



2.5. Constraints on regional averaging 

We are now going to study the asymptotic properties of the variables under the 
Lagrangian smoothing flows as ft —>■ oo. For this purpose we have to go first into the 
constraint equations in order to establish a link between the actual and the regionally 
smoothed-out initial data. 
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The constraint equations 

If we are willing to assume that Einstein's equations hold on the regions 
Bp(pj] ro) where we are averaging, then besides the integral constraint of regional mass 
preservation, we must also require that our regional averaging procedure is such as to 
respect the constraints Eqs. (|]), at least on the given scale r . Thus, we restrict the 
class of possible deformation flows to act within the solution space of the constraints. In 
other words, for each /3 the smoothed geometry should be a candidate for an admissible 
initial data set of Einstein's equations in Bp{jpy, r ). 

The divergence constraint 

Let us start by discussing the divergence constraint, which, in terms of the trace- 
free distortion tensor V ab reads: 



V 6 



^V a K = 8nGJ a , (94) 



or more explicitly 

^ [V b (V ab + £ $g ab ) - (V ab + £ $g ab ) V b In N] - 2 - V a K = 8nG J a . (95) 

If we assume that such a constraint holds in the regions Bp(pj]ro) also for the (3- 
deformed distortion tensor V ab (j3) we formally get: 

[V & (Vl b + £^g ab ) - (V* ab + £^g ab ) V b ln 

~V a K{(3) = 8nGJ a ((3) , (96) 

where the notation [....]p is a shorthand for the explicit /3-dependence of all the quantities 
within the brackets. 

According to the choice ( |9~1~D for the /^-dependent shift vector field N a ([3), (p6|) 
reduces to 

[V 6 (V* ab + £ d9ab ) - (V* ab + £,g ab ) V»lnN] p 

- 2 -V a K{f3) = 8<KGJ a {(3) , (97) 

where the matter current J a (P) is given by (p3|). 

Assuming that the smoothing flow (3 — > {g ab {f3) ) V* ab {(3),v a {l3)) is global and that 
9ab(P) g a b * s an isolated constant curvature metric, then T>* b (/3) — > 0, V b lniV — > , 
and J a (P) — > 0, uniformly, and Q9"7| ) reduces to 

V b (£vg ab )-^d a K = , (98) 

where A is the j3 — > oo limiting value of the given quantity A. Thus, the possible 
anisotropies in the gradient of the rate of volume expansion K, as seen by the observers 
executing the smoothing process in Bp(pj-,r ), are only due to the diff-warp generated 
by the Ricci flow. 
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The Hamiltonian constraint 

We are now in position of discussing the Hamiltonian constraint. By adopting 
the same regional philosophy applied to the divergence constraint, we assume that it 
holds in Bp(pj;r ) for the /^-dependent data on (E^, g a b(fl))- Explicitly, in terms of the 
trace-free distortion tensor V ab , we get: 

~wwf b{(3)9Cd{(3) ^ + £ ff^p + £ 

= l6nGp(/3) + 2A . (99) 
According to (f?H) this reduces to 

n{(3) + 2 -K\(3) - ^jf b ^)9 cd ^) IKc + £v9ac) p [D* bd + £v9u] p 

= l6nGp((3) + 2A, (100) 
which, upon taking the volume-average over Bp{j>j\ r ), yields: 

1 / g ab (P)g cd (P) IK C + £wac] p [V* bd + £,g bd ] p 



4 \ N*{(3) 
= 16ttG (p(f3)) Bg + 2A . (101) 
According to (^) we have: 

(p(/?))^ = (p), o e / ' f t WA)> ^- WA)>s > A • (102) 

Since, in the f3 — > oo limit T>l d (f3) — > 0, N(f3) — > 1 uniformly (again by assuming that 
the Ricci flow metric g a b{(3) — > g~ ab is an isolated constant curvature metric), we can 
write (|101|) as: 

= 16.G(p) Bo e / "l <K( ' 5 H-« K <™=>« + I( r y^ s ^£ j M ) if+ 2A . (103) 
Observe that 

\{T b g cd £vg ac £^g bd )^ = 2{a% , (104) 

where a 2 = ^a ab a ab is the squared norm of the shear tensor. Moreover, since g ab is 
a constant curvature metric, TZ is a constant that, in order to emphasize the regional 
nature of the averaging, we denote by IZg. Thus, we finally get: 

\(K 2 )_= 16 7 rG(p) Bo e^N /3)> ^- < ^ )>s J d/3 + 2<a 2 > F + 2A-% . (105) 
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2. 6. Summary of key-results 

Before we discuss cosmological implications, we here summarize the key-results of the 
previous sections that are relevant for applications. 

In Section 1 we have introduced the concept of a (position-dependent) system of 
geodesic ball-domains on which volume-averages of scalar functions are evaluated. We 
have explicitly shown how the averaged density of matter and the geometry of balls 
change under variation of scale. Complementary to this (Eulerian) averaging under 
variation of the ball radius, we have devised a (Lagrangian) smoothing flow that provides 
a conceptually equivalent averaging procedure. Here, we gave substance to the choice of 
the Ricci flow as a natural candidate for the smoothing of the metric. The key-results 
of this section were: 

• For scalar functions on a Riemannian 3-manifold the operations spatial averaging 
and rescaling the domain of averaging do not commute. In particular, this 
shows that averaging of (scalar) inhomogeneous fields implies the necessity of 
simultaneously rescaling the (tensorial) geometry. 

• The metric in the neighborhood of the domain of averaging is forced to rescale in 
the direction of its Ricci tensor. 

Important equations were Eq. ([44]) (Proposition 1) and Eq. (|33| ) in the context of 
(Eulerian) averaging under variation of the ball radius, and the corresponding equations 
in the context of (Lagrangian) smoothing, Eq. fl47D and Eq. fl48|) . 

In Section 2 we have put into practice the smoothing program in terms of a globally 
volume-preserving Ricci deformation flow. We stressed that the choice of this global 
normalization is technical, the key-results concerning regional averages do not depend 
on this choice. We devised a corresponding deformation flow for the material mass 
under the assumption that the total material mass within the domain of averaging is 
preserved during the deformation. We then showed that the stability properties of the 
Ricci flow entail a unique choice of the smoothing flow for the second fundamental 
form. We implemented the requirement of the preservation of the constraints under 
such deformations, and determined the optimal choice for the reference frame in which 
fundamental observers execute the smoothing procedure. The key-results of this section 
were: 

• The Ricci flow for the metric and the corresponding material mass flow link the 
global and regional averages for the material mass density and the scalar curvature. 

• The second fundamental form must necessarily be deformed according to the 
linearized Ricci flow; there are no other consistent choices. 

• The optimal choice of the smoothing reference frame is determined by compensating 
for the fictious motion of fundamental observers that is induced by the geometrical 
deformation. 
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To get a consistent picture, the constraint equations were required to be preserved 
during the smoothing procedure, i.e., the smoothing flows were required to act 
within the solution space of Einstein's equations. 



Important equations were Eq. ([52]), Eq. fl53|) and Eq. fl5"8|) for the Ricci- and 
material mass flows, and Eq. ([J7]) together with Proposition 2 for the linearized Ricci 
flow in connection with the smoothing of the second fundamental form. Instead of 
smoothing the second fundamental form, smoothing of the trace-free distortion tensor 
was suggested, Eq. (|77|) . Eq. (|8lf ) and Eq. (|9ll) determine the smoothing of the lapse 
function and shift vector, respectively. As for the constraints we note that, while the 
momentum constraint was of technical importance concerning the consistency of our 
choice of smoothing the shift vector field, Eq. fl97f ), (remember that the shift vector 
field cannot be disentangled from the averaging of the matter current density, and the 
averaging of the distortion tensor T> a b(f3) generates, according to Lott's factorization 
theorem, a geometrically induced shift), the Hamiltonian constraint is essential for the 
following applications, so we especially point out Eq. ( |101|) and Eq. (|105|) . 



In the next section we are going to discuss average characteristics in the smoothed- 
out region. In particular, we shall discuss the effect of averaging and scaling, respectively, 
on the parameters of an averaged inhomogeneous cosmology with the key result: 

• Cosmological parameters, as they are interpreted in a smoothed cosmology, are 
'dressed' by the removed geometrical inhomogeneities. 



Important equations will be Eq. ( |118| ) together with Eq. ( |1 19[ ) featuring a novel 



"curvature backreaction" effect; Eq. ( 120 ) as compared to Eq. for generalizations 

of Friedmann's equation in the actual cosmological model and the smoothed model, 
respectively. The Hamiltonian constraint may be cast into a constraint equation for 



regional cosmological parameters in the actual model (Eq. (|125|) ) and the smoothed 



model (Eq. (|129|) ), which lie at the basis of our interpretation concerning 'bare' and 
'dressed' cosmological parameters. 
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3. Cosmological Implications 



We are now going to study Eq. ( |105| ) in a cosmological setting, addressing especially 
the role of the cosmological parameters. Before we discuss them, let us first study and 
estimate the fluctuations in the rate of expansion and the scalar curvature. 



3.1. Estimating fluctuations 



In order to write ( |105|) in a form suggesting a generalized Friedmann equation, let us 
first introduce the regional variance of the distribution of K in B{jpj] ro) measuring the 
spatial fluctuations in the rate of expansion: 



K 



(K - {K) v ) 2 )_ = {k')_- {K)l 



(106) 



The second step is to exploit ( |9"8"D in order to get an estimate for &L in B(pj] r ). Note 
that the following estimates are just indicative for the sake of understanding; we shall 
comment on strategies for full estimates later. 

Recall that for any (C 2 ) function on B(pj] r ) we can write (as long as lZ(pj)rl <C 1): 



V E rl 



f dm 



/fe) + f 



+ o(rg) 



(107) 



A/to, 

'B(j>f,r ) 

where Ve is the Euclidean volume of the unit ball in M 3 , and A is the Euclidean 
Laplacian. Upon applying such a normal coordinates estimate to / = K and f = K, 
we obtain for the fluctuations of the rate of expansion: 



K 



10 



2 (6 ab d a Kd b K) ( Pj ) + ^K\p 3 ) 



+ o(r 2 ) 



which, upon substituting the expression (B3) for d a K, yields: 



^=10 



~ {S ab V c a; lac v d w; lbd ) fa) + ^k\ Pj : 



+ o{rl) 



(108) 



(109) 



where a,. 



(P) 



Inserting 

X[W)) 

B e 



£vg ac is the residual longitudinal shear, 
into the integral of (^) we also get: 



(p) 



Bo 



v B ,(P 



oo 



In) Y&L 



(P)l 



A normal coordinates estimate as in (|107|) (set / = 1) explicitly provides: 

1 



(p)bo = (p)b- 



±V-r 2 
30 '%'o 



+ o{rt) 



(110) 



(111) 



revealing the actual dependence of {p) B on the local curvature l^iPj) with respect to 
the regional average curvature IZ-g. In this connection, it is also worthwhile to discuss 
how the regional curvature in the smoothed region IZg is related to the actual average 
curvature (TZ) Bq . 
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Let us start by noting that the (normalized) (3— evolution of the scalar curvature 
obeys the following equation |T|] 



cVZ(p) 
dp 



A m TZ((3) + 2W(P) - -K(P) (Tl{(5))^ 



(112) 



where we have defined 

W(/3) := g ab g cd n ac (P)n bd (P) 



= g ab g cd [n ac + -g ac n\[TZ bd + -g bd n} 
A straightforward calculation then provides: 



n ab n ab + l -n 2 



(113) 



V a K{(3) da a q{0) + 2(K ab (P)K ab (P)) B 



-\{(n{p) - (np))B P ) 2 )B P + \w)) B0 (w)) B0 - w)> s j 



(114) 

On the closed manifold E, i.e., taking Bp — > Hp, the first term on the r.h.s. of 
Eq. (|114j) reduces to a flux through an empty boundary and the last term, specifying 
the deviation from the regional to the global average curvature, vanishes, and Eq. ( |114| ) 
may be integrated to yield: 

1 ' - <*& 



J /3=0 



2(n ab n ab 



dp 



(115) 



which shows that metrical anisotropies (i.e. lZ ab ^ 0) tend to generate a "Friedmannian" 
curvature TZ that is larger than the actual averaged curvature, whereas large fluctuations 
tend to an underestimation of 1Z with respect the real distribution of curvature. 
This curvature difference is described by a term that reminds us of the "kinematical 
backreaction" (having a similar form built from the extrinsic curvature) often discussed 
in relation to cosmological averaging, e.g., in: |§, f|. We therefore introduce as a 
notional shorthand the global curvature backreaction: 

\{(K{P) - (R(P))^) 2 K ~ m ab (P)n ab (P)h} dp . (116) 



Q 



n 



Notice that this term vanishes for a FLRW space section, so that globally Eq. ( |115| ) 
compares two constant curvature models. These models, in general, also differ by their 
global volume, which is not manifest in this equation, because we normalized the Ricci 
flow to be globally volume-preserving. 

On a regional domain of averaging, on which we concentrate in this paper, we can 
go one step further and use Eq. (|57D to express the difference between the global and 
the regional average curvature. We can then cast Eq. ( |114j) into the form: 



V, 



-2/3 



d_ 

dp 



(nP))B,v^(p) 



Vb h 



V a TZ{P) dol 



+ 



2{n a \p)n ab {p)) B0 



\{{np)-(np))B ) 2 )B, 



(117) 



In what follows we neglect the flux of curvature through the boundary of the ball, since 
we think that this term will not be of any observative relevance, at least on sufficiently 
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large portions of the Universe. A formal integration of Eq. ( |117| ) then provides the 
desired relation between the (constant) regional curvature in the smoothed model and 
the actual regional average curvature: 



2/3 



where we have introduced the regional curvature backreaction: 



Q 



R 
B 



oo 



VbM 



dp. (119) 



The integral Eq. (|118|) has the merit that it provides a transparent separation of the 
relevant terms: first, the volume effect, which is expected in this form by comparing 
two constant curvature space sections with the same matter content for which the 
regional curvature backreaction vanishes (remember that a constant curvature space 
is proportional to the inverse square of the radius of curvature, hence the volume- 
exponent 2/3); second, the curvature backreaction effect itself, which consists of a bulk 
contribution and a flux contribution through the boundary (that we neglected). Both 
encode the deviations of the scalar curvature from a constant curvature model, e.g., a 
FLRW space section. 

We are now going to relate our findings to suitable cosmological parameters by 
moving to a notation that is familiar to cosmology and accessible to the interpretation 
of observations. 



3.2. The generalized Friedmann equation 

In standard cosmology we are used to discuss cosmological parameters that are defined 
on the basis of a homogeneous-isotropic solution of Einstein's or Newton's equations 
for a self-gravitating distribution of matter. A refinement of the standard model has 
been suggested recently (0 in Newtonian cosmology, and 0, in general relativity), 
where the (global) homogeneous values of the relevant variables were replaced by their 
(regional) spatial volume-averages. For example, an averaged dust matter model in 
relativistic cosmology was found to obey a set of generalized Friedmann equations 0, 
from which we only need the averaged Hamiltonian constraint here (see also fL6f): 

6H 2 Bo - 16nG(g) Bo - 2A + (K) Bo = -Q Bo , (120) 

where we have defined, on the averaging domain Eq, the regional Hubble parameter as 
1/3 of the spatially averaged rate of expansion 9 := —K: 

Hb := if*)*, • (121) 

This form of the volume-averaged Hamiltonian constraint has the merit to isolate an 
explicit source term (the kinematical backreaction), which quantifies the deviations of 
the average model from the standard FLRW model equation. It is composed of positive- 
definite fluctuation terms 

Qf := 2(II) Bo - 2 -(lf Bo = jj<(0- (*)«,)% ~ 2(c- 2 )b ■ (122) 
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Here, I \= K = -6 and II := \[K 2 - K ab K ab ) = \6 2 - a 2 denote two of the three 
principal scalar invariants of the extrinsic curvature; the latter equality features the 



corresponding kinematical invariants expansion rate 6, and rate of shear a := y \o~ a bO~ ah 
(for irrotational flows). 

In contrast to the standard FLRW cosmological parameters there are four players. 
In the former there is by definition no kinematical backreaction, Q^ Q = 0. 

Regional cosmological parameters 

Furthermore, in the general model, we may define regional cosmological parameters as 
the following (scale-dependent) functionals B: 



and, in addition to the standard parameters (|123|) , 



n ff : = -3- ■ (124) 
These "parameters" obey by construction: 

«Sb + Q i + ^1 + V§« = 1 , (125) 

and they would all become (3— dependent functions under the smoothing flow. Eq. ( |125| ) 
furnishes a way of writing the volume-averaged Hamiltonian constraint that is best 
accessible to observational interpretations. 

However, unlike in Newtonian cosmology, where the corresponding equations have 
(apart from the definitions of the curvature and backreaction parameters) a similar form, 
it is not straightforward to compare the above relativistic average model parameters to 
observational parameters. The reason is that the volume-averages contain information 
on the actually present inhomogeneities in the geometry within the averaging domain. 
In contrast, the "observer's Universe" is a Euclidean or constant curvature model|. 
Consequently, the interpretation of observations within the set of the standard model 
parameters neglects the geometrical inhomogeneities that (through the Riemannian 
volume-average) are hidden in the average parameters of the realistic cosmology. 

Notwithstanding, we are now in position to relate the parameters interpreted within 
the standard model to the actual parameters by studying the smoothed cosmological 
model furnished by Eq. ( |105| ), which can be written in terms of the effective quantities 
Hq, (o)g, and IZ-g as follows: 

6#J- 167rG(^) F -2A + %= -Q| , (126) 

jj The stage of interpreting observations is in many cases a Newtonian cosmology. In standard 
cosmology it is common practice to introduce a frame that is comoving with a global Hubble flow; in that 
frame the constant curvature is a parameter in the background FLRW model and the inhomogeneities 
are studied within a Euclidean space section. 



Bo 
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where we have defined the residual kinematical backreaction (after smoothing) by: 

Sf:=|^-2<a%. (127) 

Thus, a "Friedmannian bias" in modelling the real (observed) region of the Universe 
with a smooth matter distribution evolving in a homogeneous and isotropic geometry, 
inevitably 'dresses' the matter density (g)g, the Hubble parameter H-g, and the 
curvature IZg with correction factors, even if the kinematical backreaction effect is 
respected. (Note that the latter is expected on a regional domain due to cosmic variance 
of the variables; see our discussion below). 

Correspondingly, an observer with a "Friedmannian bias" would interprete his 
measurements in terms of the 'dressed' cosmo logical parameters: 

tjM MGM B -a _ A r _ Kz qk _ q| 

oVqHq "-"S Oiig 

which again, by construction, obey: 

n|f + + n§ + n| K = i . (129) 

Our subsequent analysis will be focussed on discussing the actual relevance of the 
geometrical correction terms. 

The relation between 'bare' and 'dressed' cosmological parameters 

Following from our previous analysis, especially of the expansion and curvature 
fluctuations, we can collect the formulae in order to relate the 'bare' and 'dressed' 
cosmological parameters. 

Defining the fraction between the volume of the smoothed constant curvature region 
and the volume of the original bumby region, as well as the fraction of the corresponding 
Hubble parameters, by 

; «:=f| , (130) 

V Bo H B 

we can formally writefjj]: 

^Bo ^B ^B (fyBo 2/3/1 i ^ 2fo /iqi% 

where in the last equation for the fraction of the curvature parameters we introduced 
the dimensionless regional curvature backreaction parameter \x := Q§ /T^b- 

This set of equations furnishes a formal basis within which the results of this paper 
can be interpreted with respect to their relevance for observational cosmology. In the 
following discussion we comment on possible strategies for a quantitative analysis of the 
results. 

ff The denominators have to be nonzero; degenerate cases must be treated differently. Note that, e.g., 
the regional average curvature is in generic situations nonzero. 
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3.3. Discussion 

The above listed relations appear to provide a formal recipee to apply the results of 
this paper. However, it is clear that a quantitative estimate of the relations between 
'bare' and 'dressed' cosmological parameters must be based on dynamical models and 
cannot follow solely from geometrical consequences of the smoothing procedure on a 
given spatial hypersurface. 

Dynamical estimates can be subtle in a relativistic setting, since realistic models 
for the evolution of structure are well-implemented only in the Newtonian framework. 
Recently, progress has been made in estimating the kinematical backreaction parameter 
fig* in Newtonian cosmology. Putting the key-results obtained by a realistic Newtonian 
model for the evolution of structures into perspective (fflQl, f34|), it was found that, 
e.g., on a sufficiently large expanding region (of the order of several hundreds of 
Megaparsecs) , the kinematical backreaction parameter is quantitatively small, which is 



in conformity with other (including relativistic) estimates ([p5fl, ||, [[J, pEOfl , [p6fl , p9| ). 
The surprising result, however, was that backreaction can have a large influence on the 
other (standard) cosmological parameters during the dynamical evolution. Although 
the Newtonian model requires that the backreaction vanishes on the global boundary, 
we may argue in the context of the present work that, on a given space section at a fixed 
time of observation, the kinematical backreaction is quantitatively less important than 
the standard parameters, and reflects cosmic variance of the measured variables, the 
presence of which is expected on a regional domain of the Universe. Especially the work 
T~0|1 has clearly shown, however, that the values of these parameters are not related to 



their initial values evolved by a FLRW cosmology. 

According to our analysis of the effect of smoothing, we found that a similar term 
can be identified: the curvature backreaction, which we expect to play an analogous 
role as the kinematical backreaction. In line of these thoughts we therefore suggest 
to concentrate further quantitative investigations on the following, albeit at this level 
formal, considerations. 

Observe that together with the two generalized Friedmann equations in the form 
( |125| ) and ( |129| ) for the average ('bare') observables in the real manifold, and for the 
'dressed' observables in the smooth constant curvature model, 

+ + n| + uf = i , (132) 

we may consider fractions of various cosmological parameters in order to eliminate, say 
the fraction of the Hubble parameters a, and conclude on the values of the others. 
Our goal is to relate observationally determined values of the 'dressed' parameters (i.e., 
as interpreted with a "Friedmannian bias") to the actual parameters of the average 
cosmological model. 

According to the above discussion of the effect of the kinematical backreaction, 
we are encouraged to consider a subcase in which the backreaction parameters are 
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quantitatively negligible. (This is, e.g., also true when the positive-definite fluctuation 
terms compensate each other: for the kinematical backreaction we would consider 
the compensation of shear fluctuations and expansion fluctuations; for the curvature 
backreaction we would think of a compensation of fluctuations of metrical anisotropies 
with curvature amplitude fluctuations.) It should be emphasized that the following is 
an illustration and does not replace the necessity of a fully dynamical investigation. 
Approximating the 'bare' and 'dressed' kinematical backreaction parameters by zero, 
we have: 

Let us furthermore consider a region of the Universe on scales of the order of 1 
Gpc, where (also in the same spirit possibly biassed) observations of the first doppler 
peak in the CMB fluctuations at the "Friedmannian scale" ~ lOOMpc/h favour an 
approximately vanishing average curvature TZq ^ 0. If, again for simplicity, we 
approximate also the curvature backreaction parameter by zero, \i ~ 0, in the sense 
that there are curvature fluctuations present, but the two positive-definite parts in the 
backreaction term compensate each other, we would have an approximately vanishing 
average curvature also in the actual cosmological model. Then, the standard argument 
requires compensation of the actually observed matter content (including dark baryonic 
and possibly dark nonbaryonic matter components), obeying the commonly agreed 
upper bound Q-g < 0.3 with a cosmological term Q-g « 0.7. For the 'bare' parameters 
we then obtain ^b /^Bo ~ oJ u i wn i cn yields the estimate: 

«X - vfk; ; n i w 1 - < • (134) 

This (certainly oversimplified) example shows that, instead of postulating the presence 
of a large cosmological term, the 'bare' mass parameter could still acquire values close 
to one, if 'undressed', and if the volume fraction v is substantially greater than 1. The 
second relation in Eq. (|131 ) then shows, that the actual Hubble-parameter would be 



larger than the 'dressed' Hubble parameter. The estimation of v itself is intimately 
connected with the estimation of curvature fluctuations, and a quantitative (scale- 
dependent) statement about this effect is beyond the scope of this paper: it requires 
dynamical considerations. 

A quantitative estimate that gives us an idea of the order of magnitude of such an 
effect has been worked out by Hellaby |30[| . Using spherically symmetric solutions he 



compares, on some spatial scale, the volume of a FLRW space section with that given 
by a Lemaitre-Tolman-Bondi solution (see also the alternative comparison of volumes 
suggested in [jHfl ). Contrary to our approach (that compares the two models at equal 
material mass), Hellaby uses "volume matching" as proposed by Ellis and Stoeger p4[ 



(requiring the volumes to be equal). However, we expect the estimated effect to be of 
the same order which, for two realistic density profiles of a typical cluster of galaxies 
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(an over-density) and a typical 'void' (an under-density) , yields errors in the range 
10 — 30 % for these single objects, if, e.g., the spatial averages of the density profiles are 
compared with the corresponding (fitted) FLRW parameters. 
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Remarks 

Remark 1: In order to put such qualitative remarks on a firmer ground, we need to 
quantify in which sense the domains B(p; 7") peS do not differ too much from the standard 
Euclidean ball B E (0; r). To this end it is useful to introduce the Reifenberg norm of 
(J2,g a b) at the scale r, defined according to 

\\(X,g ab )\\ r ■= ~ max d G - H [B(p;r),B E (P;r)] , (135) 

where dc-H [B(p; r), B E (0; r)] denotes the Gromov-Hausdorff distance between B(p; r), 
and B E (0;r) , (recall that da-H [B(p; r), B E (0; r)] < e if B(p;r), and B E (0;r) can be 
isometrically embedded in some metric space where they have Hausdorff distance < e) . 
It is immediate to check that ||(E, gab) \\ r — > as the scale r — > 0, (simply because for 
Riemannian manifolds, the geometry of the geodesic balls approaches the Euclidean 
geometry as r — > ). Obviously, we are interested in estimating ||(£, «7 a &)|| r for finite r. 
As a practical matter, this appears to be a difficult task, owing to the rather abstract 
definition of ||(E, # a 6)|| r - However, if we assume that the scalar curvature is bounded 
below, i.e., 

M{K ab u a u b : g ab u a u b = 1} > C > -oo , (136) 

u 

then according to a remarkable result of Colding the Reifenberg norm ||(E, <7o&)|| r 
is small, for finite r, if and only if the volume of {B(p\ r)} x6 s is close to the volume 
of the Euclidean ball B E (0;r). This result provides a very useful practical criterion for 
controlling a scale-dependent averaging procedure on (E,g a b), also because once it is 
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established that ||(E, <? a &)|| ro is small, say by comparing V (B(p;r )) with V (B E (0; r )) 
at a given observative scale r , then the smallness of || (E, g^) \\ r holds at all scales r < r . 
Thus, if the Reifenberg norm ||(E,g a ft)|| r is small, i.e., if the geometry of {B(p;r)} p( zs 
does not vary too wildly with p, then the system of regional averages {(^) B ( p . r )} p es 
appears as a reliable smoothing, at every scale r, of the given matter distribution g. 

Remark 2: Even if polar geodesic coordinates suggest themselves as the most natural 
labels for the points of B(p;r), they suffer from the basic drawback that their domain 
of definition cannot be a priori estimated and strongly depends on the local geometry of 
(E, g a f,). In this connection, a much better control on the geometry of the balls B(p; r), 
and hence on {g) B ^ p . r y can be achieved by labelling the points exp~ 1 (g) e B E (0,r) with 
harmonic coordinates, i.e., a coordinate system {z 1 } such that the coordinate functions 
z % are harmonic functions with respect to the Laplacian on (E,g ab ). We can do this by 
starting from the given (Cartesian) normal coordinates y\ and look for a diffeomorphism 
on a sufficiently small Euclidean ball B E (0,r) C M 3 , 



The standard theory of elliptic partial differential equations implies that the harmonic 
functions so characterized do form a coordinate system in Be(0, r). The important 
observation is that such harmonic coordinates can be introduced on balls of an a priori 
size as soon as the manifold (E,g afe ) has bounded sectional curvature and its injectivity 
radius is bounded below. 
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